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The dual superconductivity of the vacuum in SU{3) gauge theory is investigated by constructing a disorder parameter 
which signals monopole condensation in various abelian projections and by studying numerically on the lattice its behaviour 
at finite temperature. We find that the vacuum is a dual superconductor with respect to each f/(f ) of the residual gauge group 
after abelian projection independently of the abelian projection chosen. Like in the SU{2) case (discussed in a companion 
paper) a finite size scaling analysis enables us to extract the indices of the phase transition and our analysis is consistent with 
independent determinations. 



PACS numbers: ll.15.Ha, 12.38.Aw, 14.80.Hv, 64.60. Cn 



I. INTRODUCTION 

In a companion paper |l|], which we will quote as I, 
we have presented the basic ideas about confinement 
and dual superconductivity of the ground state of gauge 
theories, and how they can be tested in SU{2) gauge 
theory. Monopoles exist in gauge theories, carrying a 
conserved magnetic charge. We have defined a disor- 
der parameter (/x) detecting dual superconductivity as 
spontaneous breaking of the U{1) symmetry related to 
magnetic charge conservation, (/i) ^ signals that the 
ground state is a superposition of states with different 
magnetic charge, a phenomenon which is denoted as con- 
densation and which implies dual superconductivity un- 
der very general assumptions. 

In SU{2) a monopole species can be associated to any 
operator in the adjoint representation, with a correspond- 
ing magnetic U{1) symmetry. Condensation can be nu- 
merically investigated for different monopole species, in 
connection with confinement, by lattice simulation at fi- 
nite temperature. The main results of this investigation 
for SU{2) were the following: 

• monopoles defined by different abelian projections 
do condense in the confined phase, or (/i) ^ 0; 

• at deconfinement (/i) — > 0; 

• a finite size scaling analysis allows to determine the 
critical index v of the correlation length, the crit- 
ical /3 and the index S by which (/i) -^ 0. The 
determination of v agrees with the ones done by 
other methods, and indicates a 2nd order phase 
transition. Also Pc coincides within errors with 
the known values; 

• all the monopole species considered have a similar 
behaviour, and show dual superconductivity. 



Our conclusion was that confinement is an order-disorder 
transition. The symmetry which characterizes the dual 
order is not fully understood, but for sure the different 
(/i)'s are good disorder parameters. 

In this paper we shall extend the analysis to SU{3) 
gauge group. The essentials are not changed with re- 
spect to SU{2). Some formal complications come from 
the coexistence of two monopole charges for each abelian 
projection (Sect. y). We have performed a systematic 
numerical investigation, for different abelian projections. 
Also for SU{3) we find dual superconductivity in all the 
abelian projections that we have considered, again indi- 
cating that the guess of ref. M that all monopoles are 
physically equivalent is correct, {/i) looks, within errors, 
the same for the two independent monopole charges of 
a given abelian projection. A finite size scaling analysis 
shows that the transition is first order. Numerical details 



and results are given in Sect. [II. 

Sect, tt^ contains some concluding remarks. 



II. THE ABELIAN PROJECTION. CONSERVED 
MONOPOLE CHARGES 

In analogy with the SU{2) case we shall denote by 

8 

cb{x) =Y,^\x)F^ (1) 

the generic local operator in the adjoint representation. 
F* = AV2, with A' the Gell-Mann matrices. We shall 
assume (f> hermitian, or 0* real in any configuration. 
It will be convenient to use the notation 

4>{x) = $ix) ■ F (2) 

for cq. (m) and for any two operators 4>i , (f>2 

8 

2tr(/.i02 = (?i-02=^0102 ■ (3) 



Any 4'{x) can be diagonalized by a unitary transfor- 
mation U{x) 

Uix)<l)ix)U^{x)^Mx) ■ (4) 

In the usual representation of the A matrices 

q^D ^ ^IF^ + ^'dF' , (5) 

where F'^ and F^ are independent hnear combinations of 
F^ and F^. We shall choose 



F" ^F^ , F'' = —F^ + — 



for reasons which will be clear below. 
We now define 



rix) = (Uix))^ F'^Uix) , 
f{x)^{U{x))^F'U{x) . 



(6) 



(7) 



U{x) is defined as the matrix which diagonalizes 4>{x). To 
eliminate ambiguities the eigenvalues can be ordered in 
decreasing order. U{x) is determined up to an arbitrary 
matrix Ud{x) on the left U{x) ~ Ud{x)U{x), with Ud = 
exp (iaF^ + if3F^^, i.e. up to a residual U{1)'^. From 
eq. (P) and ffl), in the usual representation of Gell-Mann 
matrices 



1 / ^D + 2^D 



»_D 



2a/3 





^d-V>''d I , (8) 

-2^'h - f'h 



with ifl) >0, ip''jj> 0. 

The gauge transform U{x) is singular at the sites where 
either (/jJj = 0, and 



or cpf) = 0, and 



1 
2\/3 I -2 



2 

^D I -1 

0-1 



(9) 
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In both cases (J)d has two equal eigenvalues. 
The two field tensors 

F^l;.' = ^tr l^r'^ ■ G,. - '-r''- 



(11) 



[D^.r\ D,r'']\ 



are the analogous of the 't Hooft's tensor Q in SU{2). 
Like in SU{2) the bilinear terms in A^A^ cancel. In the 
abchan projected gauge /"(s) = i^° and f^{x) = F^ are 
x independent, apart from singularities, and therefore in 
the domain in which U{x) is regular 



A"' 



d^A^^' 



d,A 



a.b 



(12) 



The cancellation of the bilinear term A^^Ay between 
the two terms of eq. ( pT| ) is not automatic in SU{i) for 
arbitrary choice of /", /'', as it was in SU{2), and only 
works if /", /'' belong to f/(l) in the breaking SU{3) -> 
SU{2) X [/(I), which is the case for the choice of (g). Also 

the choice ^F^ — ^F® for /" or /^ would be legitimate. 

As in SU{2) F^^;}'*, the dual tensor to F^;j>, define 
two magnetic currents d'^F'^^i^* = j"'^, which are con- 
served. The theory has two conserved magnetic charges, 
M", A/^. Monopoles exist at the points where U{x) is 
singular: where ipjj = the monopole field is directed as 
F^ = diag( 1/2 , —1/2 , ), where (p'j^ is zero it is directed 
as F^" = -(1/2)F3 -I- iV3/2)F^ = diag( , 1/2 , -1/2 ). 
As in SU{2) we shall investigate the invariance of the 
ground state with respect to these magnetic [/(l)'s, in 
connection with confinement. 

On the lattice, we shall define the abelian projected 
fields as follows. In the abelian projected representation 
we write the generic link [/^(rt) in the form 



Uf.{n) 



■F^^i{VaF''+ViF'') 



(13) 



with F± a superposition of generators belonging to non 
zero eigenvalues of the Cartan algebra. 

Eq. (O) is easy to prove: it is a trivial consequence of 
the Baker- Hausdorff formula. Like for SU{2), the abelian 
part of a product is the sum of the abelian parts of the 
factors, to order a^ [a is the lattice spacing). The abelian 
magnetic fluxes through plaquettes, one for each C/(l), 
can be defined, and are identically conserved. The disor- 
der parameter is again 



(m-^) = 






(14) 



^a,6^^^^Tr{n,o(n,t)-nf/(n,i)} , (15) 



T'a,fc 



(10) Hio' ("-j^) is obtained from nio(n, i) — Ui{n,t)Uo{n + 



i, t) {Ui{n, t + 1))'' (C/o(n, t))'' by the change 



(16) 



giAf,(«,i?)/''.''(n+J,t)giA(«+i,y)/'''''(«-HJ,t) 



where 



A^'{n,y) = Af{n,zj)+VAin,y) 



(17) 



is the vector potential produced by a monopole. The 
proof that (/i"^'') creates a monopoles of the correspond- 
ing type is exactly the same as for SU{2). 

Also for SU{3) instead of (/i"''') it is convenient to de- 
termine 



p''■''^J^log{^l^■') = {S)s-{S + AS) 



S+AS 



(18) 



as a function of /3. We do that on an asymmetric lattice 
N^ X Nt {Ns 3> Nt ) which provides the static thermal 
equilibrium at T = l/a(/3). 

The deconfining transition is known and has been stud- 
ied using the Polyakov loop order parameter Q . We will 
investigate if going from deconfined to confined phase 
monopoles do condense to produce dual superconductiv- 
ity, (/i"''') will be the disorder parameters. 



more noisy. Thermalization was checked by monitoring 
the action density and the probability distribution of the 
trace of the Polyakov loop. The number of measurements 
was at least SOOrc, where tc is the correlation time of 
the considered set of data. 

For A± we use the Wu- Yang's parameterization; we 
have also checked numerically that Dirac's form gives 
similar results, as expected. 

In terms of p 



0{n)^n,j{n) 



(20) 



= U,{n)U,{n + i) {U.,{n+j))^ {Uj{n))^ ; 

O is the "butterfly" F 

0{n) = F{n) (21) 

- U^{n)Uy{n + x) {U^{n + y))^ {Uy{n))^ 



iP' 



a.b\ 



exp 




(/3')d/3' 



(22) 



III. NUMERICAL RESULTS 

We determine the temperature dependence of p on a 
lattice N^ x Nt (Ns ^ Nt), with N = i and N^ ranging 
from 12 to 32. For the reason discussed in |l|, we use pe- 
riodic boundary conditions in the spatial directions and 
C*-boundary conditions [pi in the time direction. As in 
SU{2) we diagonalize an operator O belonging to the 
group O = exp(i$"A"/2) and we identify $° by ordering 
the imaginary part of eigenvalues in decreasing order. 

As for the SU{2) case, we study the following projec- 
tions: 

• O is connected to the Polyakov line L{n,t) = 
^^irt^Uo{n,t')Ulz}oUa{n,t') in the following way: 

0{n, t) = li^L-^Uoin, t')L*{n, mV^^Uoin, t') ; (19) 

(Polyakov projection on a C*-periodic lattice); 

• O is an open plaquette, i.e. a parallel transport on 
an elementary square of the lattice 



Eq. ( p^ ) implies that if the dual C/(l) symmetry defined 
by some abelian projection and by some abelian genera- 
tor of the gauge group is related to colour confinement, 
in the thermodynamic limit Ns —^ oo the corresponding 
p stays finite in the strong coupling region (/3 < /3c) a-nd 
goes to — oo linearly with Ng in the weak coupling region 
{/3 > Pc)- In the critical region, the abrupt decline of 
(/i) is signaled by a sharp negative peak of p; the value of 
p in this region must behave as a function of Ns as pre- 
scribed by the finite size scaling theory of pseudocritical 
systems. 
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FIG. 1. p vs. /3 for different abelian projections. Lattice 
12^ X 4, abelian generator F^. 



The trace of F is the density of topological charge. 

The simulation was done on a 128-node APE Quadrics 
Machine. We use an overrelaxed heat-bath algorithm to 
compute the Wilson term of eq. (|l^), and a mixed al- 
gorithm as described in our previous paper p] for the 
monopole term (S + AS)s+as- Fax from the critical 
region at each /3 we sampled over 4000 termalized con- 
figurations, each of them taken after 4 sweeps. The er- 
rors have been obtained by using a Jackknife method on 
binned data, as discussed in W. In the critical region a 
higher statistics is required. Typically the Wilson term is 



Fig. |l| shows the typical behaviour of p for different 
abelian projections, for a lattice 12^^ x 4. As abelian gen- 
erator we used F^. The negative peak occurs at the ex- 
pected transition point, Pc @- Below Pc the different 
projections are equal within errors, suggesting that dif- 
ferent monopoles behave in the same way. 

We have investigated also whether at fixed abelian pro- 
jection the profile of p depends on the U{1) magnetic 
subgroup. Fig shows the profile of p corresponding to 
F^,F^ and F^'^ in the Polyakov projection on a 12'^ x 4 
lattice. No appreciable differences can be seen between 



different choices. This is an indication (confirmed also 
by simulations on larger lattices) that monopoles defined 
with respect to different abelian generators behave in the 
same way in the 5'f/(3) vacuum. This is also true for the 
other abelian projections we have investigated (see fig. |^). 



In the strong coupling region at low /3's_p seems to con- 
verge to a finite value (cfr. fig. ^). Eq. ( p^ ) then implies 
that (/i) 7^ in the infinite volume limit in the confined 
phase for these /3's. Hence monopoles do condense in this 
phase. 
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FIG. 2. p vs. f3 for different abelian generators. Lattice 
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12 X 4, Polyakov projection. 



FIG. 4. p vs. /3 in the strong coupling region for lattice 
sizes Ni, X 4. Polyakov projection, abelian generator F"'' . 
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FIG. 3. /9 as a function of /3 for different spatial sizes at 
fixed Nt = 4. Polyakov projection, abelian generator F'' . 



FIG. 5. p vs. Ns {Nt = 4) at /? = oo in the Polyakov 
projection with abelian generator F^. Data are obtained by 
numerical minimization of {S + AS)s+as- 



Since different abelian projections and different abelian 
generators give indistinguishable results, for the sake of 
simplicity we shall only display the Polyakov projection 
and the abelian generator F^ in the following figures. 

Fig. H shows the dependence of p on N^. The qual- 
itative behaviour does not change when we enlarge the 
lattice size. 

We now analyze the dependence on Ns in more detail. 



In the weak coupling region, we can evaluate p pertur- 
batively. The path integral is then dominated by the clas- 
sical solutions of the equations of motion for the gauge 
variables and we have 



p — > minis'! ~ minis' -|- AS} 

/3— oo I U U 

= - minis + AS} , 

u 



(23) 



since inmu{S} = 0. 

In other systems, where the same shifting procedure 
has been appUed and studied, this asymptotic value has 
been analyticaUy calculated in perturbation theory with 
the result 



P 



-cN, 



(24) 



where c and d are constants, i.e p goes linearly with the 
spatial dimension. 

In SU{3) we are unable to perform the same calcula- 
tion and we have evaluated the minimum min[/{S'-|- AS"} 
numerically. Detail about the followed procedure have 
been discussed in [Q . Here we note that due to the single 
precision of the APE Quadrics Machine, the estimation 
of the minimum of 5 + AS* for the biggest lattice is more 
noisy than in the SU{2) case. 

The result is shown in fig. S for the Polyakov projec- 
tion. It is consistent with the linear dependence of eq. 
(p3) with c ~ 2 and d ~ —12. Thus in the weak coupling 
region in the thermodynamic limit p goes to — oo linearly 
with the spatial lattice size and 



(m> 



N^ 



A^{-cN^+d)P ^ Q^ /3>/?c 



(25) 



The magnetic U{1) symmetry is indeed restored in the 
deconfined phase. 

The behaviour of p in the critical region can be inves- 
tigated by using finite size scaling techniques. We know 
that the transition is weak first order with a behaviour 
which is difficult to distinguish from that of a second or- 
der transition. 

By dimensional argument 






(26) 



where a and ^ are respectively the lattice spacing and the 
correlation length of the system. 
Near the critical point, for (3 < (3c 



£,^{l3c-py 



(27) 



where v is some effective critical exponent. In the limit 
Ns 3> Nt and for a/^ ^ 1, i.e. sufficiently close to the 
critical point we obtain 



(28) 



(/.)=iV7*>$(iVy^(/3c-/3),0,0) 
or equivalently 

P^=f(Nl/^{Pc-l3)) ■ (29) 



N 



ji/^ 



The ratio p/Ng is a universal function of the scaling 
variable 



X = N]''' {Pc - P) 



(30) 



For a pseudocritical behaviour, we expect v — 1/3. Using 

also [3c{Nt == 4) == 5.6925 §], we can plot p/Nl'" as a 
function of x. 

If we perform such a plot, we find that the scaling 
relation (|2^) does not hold. Such a scaling violation is 
due to finite size effect. A relationship more appropriate 
than (H) is 

/ (n]/^ {(3c - P)) + HNs) , (31) 



where ^{Ng) parameterizes finite size effects. If we as- 
sume that these effects are not critical Q then $ is given 

by 



^{Ng 



a 



(32) 



where a is a constant. This parameterization is correct 
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FIG. 6. Quality of scaling in the Polyakov projection at 
Nt — 4. Abelian generator F^. 



Fig. ^ shows the quality of the scaling for a = 190. 
Our estimate gives v = 0.33 ± 0.07 and a = 190 ± 20. 

In the thermodynamic limit in some region of /3 < Pc 
we expect 



ifi) ex iPc - Pf 



which implies 



A^.. 



i/u 



(33) 



(34) 



Using eq. ( 34 ) it should be possible in principle to de- 
termine i>, S and Pc- Our statistic is not enough accurate 
to perform such a fit. However, we can determine S using 
as an input Pc, v, which are known, by parameterizing 
p in a wide range by the form 



TV, 



\lv 



8_ 

X 



a 



(35) 



where c is a constant, as suggested by fig. g. 

Our best fitFl^ — 0.54±0.04 for the Polyakov projection 
and compatible results for the other projections. The y^ 
is order 1. 

This concludes our argument about the thermody- 
namic limit [Ng — > oo). The deconfining phase transition 
can be seen from a dual point of view as the transition 
of the vacuum from the dual superconductivity phase to 
the dual ordinary phase. That feature seems to be in- 
dependent of the abelian projection and of the abelian 
generator chosen. 
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IV. CONCLUDING REMARKS 

Like for SU{2), also for SU{3) gauge theory we have 
found evidence that transition to deconfinement is an 
order-disorder transition, the disorder parameter being 
a condensate of magnetic charges. A finite size scaling 
analysis of the system gives critical indices compatible 
with a first order transition, in agreement with determi- 
nations done by other methods [0. 

Of course we have investigated a limited number of 
abelian projections: like in SU(2), however, the indica- 
tion is that physics is independent of that choice. 

An interesting issue would be to investigate if the 
mechanism is the same in the Nc = cx) limit. As a con- 
sequence also in the presence of dynamical quarks the 
behaviour should be similar, as well as the symmetry 
pattern and the disorder parameter should be the same. 
Investigation in this direction is on the way. 
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